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ABSTRACT 

The resolvability of point sources of incoherent light is 
analyzed by quantum detection theory in terms of two hypothesis- 
testing problems. In the first, the observer must decide whether 
there are two sources of equal radiant power at given locations, 
or whether there is only one source of twice the power located 
midway between them. In the second problem, either one, but not 
both, of two point sources is radiating, and the observer must 
decide which it is. The decisions are based on optimum processing 
of the electromagnetic field at the aperture of an optical instru- 
ment. In both problems the density operators of the field under 
the two hypotheses do not commute. The error probabilities, deter- 
mined as functions of the separation of the points and the mean 
number of received photons, characterize the ultimate resolvability 
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of the sources. 



Two point sources of light appear as one when they are very close; one 


function of a telescope or a microscope is to separate their images to the 
point of distinguishability . How well it does so is measured by its resolving 
power. When all the aberrations of the lens system have been eliminated, 
diffraction of the light at the aperture of the instrument remains to spread 
the images and cause them to overlap. According to the commonly accepted 
Rayleigh criterion, the two images are said to be resolved when the peak 
illuminance of the diffraction pattern of one falls on the first minimum of 
the diffraction pattern of the other. Alternative measures of image resolva- 
bility and instrumental resolving power have been proposed, and methods such 
as apodization for maximizing them have been studied. * These strive 
essentially to raise as far as possible the likelihood that an observer will 
see two close point sources as indeed two rather than one. However ingeniously 
his optical system may be designed, his perception is subject to error because 
of the stochastic, quantum-mechanical nature of the light. By studying reso- 
lution from the standpoint of hypothesis testing, we can determine how the 
observer can decide most reliably whether there are two sources or only one, 
and in this way we can bring out the fundamental limitations on the resolva- 
bility of two luminous points. 

Two such decision problems will be studied in this paper. In the first, 
the observer is to decide whether a single point source of known power is 
present in an object plane, or whether two sources are present, each emitting 
half the power. In the second, he is to decide which of two point sources is 
radiating during a certain interval, only one being allowed to radiate at a 
time. The sources, at known locations, emit incoherent, quasimonochromatic 
light of given spectral density. The decisions are to be based on the electro- 
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magnetic field at the aperture A of an optical instrument during a fixed 
observation interval (0, T) . The instrument is to process that field in such 
a manner that the decisions can be made most reliably. The minimum attainable 
probability of error characterizes the resolvability of the sources. 

The optimum processing of the aperture field is determined by quantum 

[3-5] 

detection theory . It requires us to find the eigenvectors and eigen- 

values of the operator pi - Apo» where Pq and pi are the quantum-mechanical 
density operators of the aperture field under the two hypotheses in question, 
and A is a constant. In neither of our problems do the density operators 
commute. Except for choices between pure states, there are few physically 
significant pairs of noncommuting density operators for which exact sets of 
eigenvalues and eigenvectors have been calculated. In order to obtain exact 
solutions, it has been necessary to assume the absence of background light, 
errors in the decisions arising only because of the quantum nature of the light 
from each source. We determine in this way an absolute limit to their resolu- 
tion. 

In attacking each problem we must first find an expansion of the aper- 
ture field in terms of which the density operators take the most convenient 
forms. Then the eigenvalue equation is written down in the coherent-state 

[g] 

representation, which leads to especially simple eigenvectors, constructed 
from the eigenvectors of p 0 and p^ individually. The vanishing of a deter- 
minant of the coefficients of these eigenvectors gives us a quadratic equation 
for the eigenvalues. From the eigenvalues and eigenvectors the probabilities 
of error can be calculated. They depend on the separation e of the two sources 
on the form and size of the aperture, and on the total average number N g of 
photons received during the observation interval. 
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1. Two Sources or One? The Aperture Modes. i . 

The observer is to choose between two hypotheses: (Hp) a single point 

source with a total radiant power P s is present at the origin u = 0 of an object 
plane at distance R, and (H^) two point sources, each with radiant power % P s , 
are present at points u = - % e and u = + h. § in that plane. The sources 
radiate quasimonochromatic incoherent light with a mean angular frequency ft , 
wavelength X = 2ttc/0 (c is the velocity of light) , and a spectral density 
X(w - ft) whose width 2nW is much less than ft. 

The decision is to be based on observations of the electromagnetic field 
at the aperture A of an optical instrument during an interval (0, T) . The 
aperture plane lies parallel to the object plane. The field can be treated 
as a scalar function ip(r, t) of aperture coordinates r and time t, provided — 
as we assume — that the sources are close enough together that the rays from 
them are paraxial, | e | /R << 1, and W << ft. The field is broken into its 
positive-frequency part ij/ + ^(r, t) and its negative-frequency part ^ (r, t) , 

<J>(r, t) = i|/ + \r, t) + <f/ ^(r, t). 

Classically, because the sources radiate natural, incoherent light, ij/ + ^(r, t) 
is a complex Gaussian spatio-temporal random process. Quantum-mechanically the 
field is described by its density operators pp and pj under the two hypotheses. 

r q 1 

In the coherent-state or P-representation these have Gaussian forms 

depending only on the mutual coherence functions of the light, which are the 

[9] 

quantum-mechanical averages 

h Tr Pi t : ) i/ +) (r 2 , t 2 ) = 

(r x , r 2 ) x(ti - t 2 ) expt-iftUi - t 2 )], i = 0, 1. (1) 
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Now t) is an operator in the Hilbert space of states of the field, 

<p( \r, t) is its Hermitian conjugate, and Tr stands for the trace; x( T )> the 
temporal coherence function, is the Fourier transform of X(w) , normalized so 
that x (0) = 1. 

The spatial coherence function <p s ^^(rj, r ) is proportional to the 

spatial Fourier transform of the radiance distribution B(u) in the object 
[9] 

plane. Under hypothesis H 0 B(u) = P s 6(u), and the spatial coherence func- 

tion 

9s°’ )( E 1 , F 2 ) = (E s /AT) Frx Fr 2 * (2) 

is constant except for the ubiquitous Fresnel factors Fri and Fr 2 *, which are 
defined as 

Fr^ = exp(ikr^ 2 /2R) , k = 2 tt / A , i = 1, 2. (3) 

Here E s = P s AT/4ttR 2 is the total average energy received at the aperture, 
whose area is A, during the observation interval. Under hypothesis Hi 

B (u) = h P s [6 (u - h e) + 6 (u + % e) ] , 

and the spatial coherence function is 

9g 1 \r 1 , r 2 ) = h (E S /AT) {exp[ike«(ri - r 2 )/2R] + 

exp[-ike'(r! - r 2 )/2R]} Fr! Fr 2 * 

= (E S /AT) cos u(xi - x 2 ) Fri Fr 2 *, 

~3 = ( x j> y j)» j = 1,2, y = ke/2R, e = |e|. (4) 

The x-axis has been set parallel to the line between the two sources. 

The optimum strategy for processing the incident light and making the 
decision is most expediently derived from a decomposition of the aperture 
field \r, t) into a countable set of spatio-temporal modes, represented 

by products of spatial mode functions hpCr) and temporal mode functions y m (t) e 


ifit 
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The y m (t) are taken as the eigenfunctions of the temporal coherence function 

V 

x(t), as given by the integral equation 

X m = T_1 I X(t x - t 2 ) Y m (t 2 ) d ^2; (5) 

Jo 

they are orthonormal with respect to the observation interval (0, T) . The 

eigenvalues Xm> which sum to 1, give the average fraction of the light going 

[9] 

into each temporal mode. When WT >> 1 there are roughly WT such modes con- 
taining a significant portion of the light. By virtue of (5) the temporal 
modes are statistically independent, and the density operators po and pj can 
be factored into products, of which each factor refers to a different temporal 
mode; we have here a spatio-temporal counterpart of the Karhunen-Loeve expan- 
sion. For simplicity we assume at first that only a single temporal mode is 

-h 

excited; its mode function can be taken as Yi(t) = T , with xi = 1« Later 
we shall without difficulty extend our results to a large number (WT >> 1) of 
statistically independent temporal modes. 

The spatial mode functions n p (r) are orthonormal over the aperture, 

hq*(r) n p (r) d 2 r = <5 pq . (6) 

[9] 

In the simplest binary detection problems these can be simultaneously 

eigenfunctions of both spatial coherence functions 9s^( r l> J 2 ) > i = 0, 1, and 
the spatial modes are statistically independent under both hypotheses Hq and 
Hj , but that is impossible here. 

In order to limit the number of spatial modes that must be considered, 
we postulate that the aperture A is symmetrical with respect to the x- and y- 
axes. Our results will be illustrated for a rectangular aperture having one 
side of length a in the x-direction and centered at the origin. The two eigen- 
functions of 9g^(ti, r 2 ) , which is proportional to cos pxj cos yx 2 + 
sin yxj sin yx 2 , are 
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ni (r) 
P2(r) 


Cj cos yx Fr, 
C2 sin yx Fr, 


( 7 ) 

( 8 ) 


where Fr is again a Fresnel Factor, and Ci and C2 are normalization constants. 
We call m (0 the cosine mode, n2( r ) the sine mode. The associated eigenvalues 
hj and h 2 , defined by^^ 


h p hp(ri) = A' 1 | r 2 ) n p (r 2 ) d 2 r 2 > 


( 9 ) 


are for the rectangular aperture 


hj = %(1 + sine a), h 2 = k(l - sine o), 

a = ya/r = ea/XR. 


( 10 ) 


They determine what fraction of the light goes into each mode under hypothesis 
H 1# and h^ + h 2 = 1 . The remaining eigenfunctions of <Pg^(ri, 52) are ortho- 
gonal to m(r) and ri2(j) and have zero eigenvalues. 

Under hypothesis H 0 all the light goes into the planar mode no'(t) = 

A Fr, but this is not orthogonal to (r) . We therefore introduce what we 
shall call the zero mode, 

no(r) = (C 0 ' + C 0 " cos yx) Fr, ( 11 ) 

choosing C 0 ' and C 0 " so that no( r ) I s both orthogonal to m(t) and normalized. 
As an even function of x it is also orthogonal to n2(j)> which is odd in x. 

The remaining spatial modes are generated by the Gram-Schmidt procedure 
in such a way as to be orthogonal to no(^)> ni(r), and n 2 (r). They are eigen- 
functions of both <p s ^(ri, 12 ) and (ji , £2) ,but with zero eigenvalues. As 

none of them is excited by the incident light under either hypothesis, we can 
disregard them henceforth. 

(+) 

The aperture field ip (r, t) is expanded in a series of spatio-temporal 
modes Hp(r) Y m (t) with expansion coefficients proportional to^^ 


7 



a pm = U^)~ h I I n p *(r) y*(t) </ +) (r, t) d 2 rdt, (12) 

where ft is Planck's constant h/27r. Quantum-mechanically a pm is the annihila- 
tion operator for photons in the (pm) mode. The density operators pq and pi 
of the field under the two hypotheses depend only on the covariance matrices 
of these coefficients; the elements of the covariance matrices can be 
expressed through (1) in terms of the mutual coherence functions of the field. 


(i) „ + 

<P = Tr p . a a„ 

pq,mn 1 qn pm 


(T/fifi) x m 

A m mn 



,Ci). 


n p *(ri) <P S (ri, r 2 ) n q (r 2 ) d^r 1 d z r 2 , (13) 


where we have used (5) . Because we are for the time being treating only a 
single temporal mode, yi(t), we put X} = 1 and drop the subscripts m and n 
referring to the temporal modes. If we now use our mode functions from (7), 
(8), and (11) and the spatial covariance functions in (2) and (A), we find 
that the covariance matrices are, with p and q equal to 0, 1, or 2, 


= N s y P y q a - «„,) a - «„,), 


pq 

,U> = 

pq 


9';; = N s h p <5 pq (1 - 6„n) 


p2' 

5 p0 ; 


q2' 


(1A) 

(15) 


where 6 pq is the Kronecker delta, equal to 1 for p = q and to zero for p f q. 


N s = E s /hf2 


(16) 


is the mean number of photons received. 



y± = A / nj[(r) d 2 r. 


and hj and h 2 are the eigenvalues defined by (9) 
ture, 

yj 2 = 2 sinc 2 (% a)/(l + sine a), 
y 0 2 = 1 - yi 2 , a = ea/XR, 


For the rectangular aper- 


(17) 
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and hj and h 2 are given in (10). The density operators in the P-representation 


are 


[ 8 ] 


with 


Pi 



(d 2 a k /-n) , 


a = (a 0 , <*!, a 2 ), a k = a kx + ia ky> 


d 2 ai 


da kx da ky ’ 


Pi (a) 


7 r 3 jdet cp (i) j 1 exp 
i = 0, 1. 



p=0 q=0 



(18) 


Each integral is taken over the entire complex a-plane. The density operators 
Pq and pi do not commute. 
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2. Two Sources or One? The Optimum Strategy. 


According to quantum detection theory the optimum strategy for choosing 

[3-5] 

between two hypotheses requires measuring the projection operator 


n - Kr^vl U(n nr ) > U9) 

n,r 

where U(x) is the unit step function, and | n n ) is an eigenvector of the 
operator pi - Apg with eigenvalue n nji . given by the operator equation 

(p i “ a po) h nr > = n nr |n nr > ; ( 2 °) 

the sum in (19) is over all the eigenstates. (The paired subscripts, which 
are non-negative integers, are used for later convenience.) The result of 
measuring J] is either the number 0, which brings the decision for hypothesis 
Hq , or the number 1, which brings the decision for Hj . Here A is a parameter 
depending on the decision criterion. If the average error probability is to 
be minimum, A is the ratio ?/(l - £) of the prior probabilities of hypotheses 
Hq and H^ . If the Neyman-Pearson criterion has been adopted, A is set to 
yield a pre-assigned false-alarm probability. The average probability of error 
is 


P e = ? Qo + (1 “ O (1 - Q d > 


= (1 - C) 


-En u 


n,r 


(n ) 

nr nr 


( 21 ) 


the false-alarm probability that Hi is chosen when Hq is true is 


Q ° = E < Vl p olV> u(n nr } ’ (22) 

n , r 

and the detection probability — the probability of saying there are two sources 
when there really are two — is 
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(23) 


Qd 


E (n |p«|n ) u(n ). 

nr' I' nr nr 


n, r 


Both Qq and are functions of A. 

Our principal task is to solve (20). By expressing the eigenvectors 
|n k > in terms of the coherent states as^^ 

F (a*) = (a | n ) exp(% j oc | 2 ) , 

nr ~ nr 

|a) = |a 0 ) |aj> ja 2 >, |a| 2 = |a 0 | 2 + laj 2 + |a 2 | 2 , (24) 

[3] 

we can write (20) as an integral equation 

/ [ R l(§*> y) " AR 0 (|3*, Y)] exp ( - | y | 2 ) F nr (y*) J"J (d 2 Y i /tr) 

J i=0 


= V F nr ( ^> 


where 


R i(ft*, y) = <§|pily> ex P 3 s( [ § 1 2 + |y| 2 ) = 

2 2 

det(I + V (l) )~ 1 exp^^ ftp* t (I + ? (i) ) _1 ? (l) ] pq Y q , 

p=0 q=0 

i = 0, 1, 

with I the 3x3 identity matrix. Here, by (14) and (15), 

R l(ft*, Y) = (1 - Vj)(l - v 2 ) exp(v 1 B 1 *Yi + v 2ft2*Y2)» 


(25) 


(26) 


v j = N j /(N j + N j = h j N s , j - 1, 2, 


(27) 


and 


R o(§*> y) = (i - v 0 ) exp[v 0 (y 0 B 0 * + yiBi*)(y 0 Yo + yiYi)J» 


v 0 = N S /(N S +1). 


(28) 


A typical eigenvalue of the density operator pj is the probability 
= (1 - vi)(l - v 2 ) V] n v 2 r 


(29) 


of finding n photons in the cosine mode t|}(r) and r photons in the sine mode 
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n 2 (r) ; the associated eigenvector is, in the coherent-state representation, 
proportional to Y* n Y* r * For r > 0 this is also an eigenvector of pj - Ap 0 , 


F nr ^Y*) = ( n! r! ) 2 Yl° Y* r » 


n = P , 
nr nr * 


r > 0, 


because p n ln ) = 0. This means that the observer should measure the number 
of photons, or the energy, in the sine mode, which can be done without dis- 
turbing the other modes; and if he finds any photons there at all, r > 0, he 
chooses hypothesis Hj . We need, therefore, to consider further only the 
strategy to be followed when no photons are found in the sine mode (r = 0) . 

The density operator p 0 has the eigenvalues 

P ^ 0) = (1 - v 0 ) v 0 n , (31) 

which are the probabilities of finding various numbers of photons in the 
planar mode no * C H") 5 the associated eigenvectors are in the coherent-state 
representation proportional to (yoYo* + yiYl*) n . As a solution of the integral 
equation (25) we therefore try 

F n0 (y *> = x 0(y<)Y0* + YlYi*) n + XiY^ n , (32) 

with x 0 and xj constants yet to be evaluated. After substituting it and 
carrying out the integration, we find by the method outlined in the appendix 

( x i + yi n x o) e* n - A(xiyi n + x 0 ) p^(yoe 0 * + yi^i *) 11 

= n n0 [x 0 (yoBo* + yi3i*) n + Xi3i n ] . (33) 

Equating the coefficients of the two types of terms, we obtain the linear homo- 
geneous equations 


<P ^0 > - n n0 ) *1 + y >" p n0 ) = °- 

- A yi n P< 0) X1 - <AP^ 0> + n n0 ) x 0 . - 0, 
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and setting the determinant of the coefficients of xq and xi equal to zero 
yields a quadratic equation for the eigenvalues n n Q, 

9 n 

(35) 


"„0 - (P nJ’ - AP n 0,) "nO ' A? n0 ^ ’ °' 


For n = 0 the two eigenvalues are zero and 
n _ p(D A p (0) 

noo ~ Poo ” ap 0 


(36) 


Whether noo is positive or negative depends on A; for A = 1, noo < Of the 

remaining pairs of eigenvalues one is positive and one is negative, n^Q^ > 0, 
n^Q^ < 0. Putting these values back into one or the other part of (34) enables 
the ratio xj/xq to be determined. The values of xq and x^ are then found from 
the normalization requirement 

nr A . 

i) 


Kol’W * 1 ■ 


fffw e »p(-hl 2 ) |~J (<i 2 Yi/»: 


n 


= n ! (xj + 2yi x 0 X! + x 0 ). 


By using (30) and 


00 oo 

(i) 

; = i, 


lulu nr 

n=0 r=0 

we can write the average error probability from (21) as 


(37) 


(38) 


= 4 Q 0 + (1 - 5 ) Q d = 


(1 - c) l min [P qo } > Ap£ U; ] + ^ (P. 

* n=l ' 

To find the false-alarm probability we use (22), in which the sum can be 

restricted to terms with r = 0 because 


,( 0 ) 




■ (1) - n (+) ) 
nO nO ' 


(39) 


(n pn n ) = 0, r > 0. 

nr u nr ’ 


( 40 ) 


We find for each term 


( "nol p «l"nO > ’ 
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jjj F nO ( $* } R ° ( §*’ ^ F nO ( ^ ) ex P<~l §1 2 -M 2 > Y\ < d2 ^i d2 YiA 2 ) 

= n!(x 0 + x 1 y 1 n ) 2 P^°\ (41) 


whereupon the false-alarm probability becomes 

( ( 0 ) < x o + x iyi n ) 2 


Qo - 


+ P$ 0) - AP$ 0) ) 


n>0 


x l 2 + ZxQXjy ! 11 + x 0 2 
- E p n 0) - ’£><1 - n 2n ))/ Rn + pS 0) U(P<;> - AP<°>>, 


n>0 


where 


,(1) 


,(0). 


,(D .(0) 


2n N A 


(42) 

(43) 


*n = ^ P n0 - A V } + 4AP nO C (1 ' )] * 

By using (39), (42), and rather much algebra, we find for the detection pro- 
bability 

Q d 


L (l) x o 2(1 yi ^ (1) (0) (1) 

P„n — Pjo UCAPr^ - P ( 00 ) 


n>0 


nO 2 o n , 2 

x o + 2x 0 x 1 y 1 + x 1 z 


= 1 -Y1 ° (1) rA ” (o) - 2n '” (1) " (+) ' 


nO t AP n yi - < P n0 ’ Ko^n 

- p5J } u(ap$ 0) - P^ } ). 


(44) 


If thermal background light is present, the covariance matrices in (14) 
and (15) have additional terms C/V'Spq, where cJf is the mean number of background 
photons per mode. The eigenfunctions F (y*) are sums of trinomials in yq*, 
Yl*> Y2* with various powers, but the equations for their coefficients and the 
determinantal equation for the eigenvalues are much more complicated. 
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Multiple Spatio-temporal Modes 


Now we suppose that the field is divided among a great many statistically 


independent temporal modes, each of which bears three spatial modes no(r), 

rji (r) , and n 2 (t). We shall show that (39), (42), and (44) still hold, provided 

that the probabilities an ^ *”n^ are replaced by Poisson probabilities. 


= N„ n exp(-N s )/n! , 


pff - rf- 


Let us consider the finite number v >> WT of temporal modes having as temporal 
eigenvalues Xi, X v , the v largest eigenvalues of x( T ) • Later we let 

v and WT go to infinity. 

The first step in the decision strategy is to determine whether there 
are any photons in any of the v sine modes Y m (t) e . If so, hypothesis 

Hj is chosen at once, for we know that this would be impossible under hypothesis 
Hq . The only eigenvectors of Pi - Ap 0 that need to be considered further, as 
before, are those in which there are no photons in any of the v sine modes. 

We add a subscript m to the complex variables y and £ to indicate the 
temporal mode they refer to, and we recognize that the coherent-state represen- 
tations Rq and R} of the density operators pq and pj are now products of functions 
like those in (27) and (28), with a factor for each temporal mode. The eigen- 
functions of the integral equation now have the form, replacing (32) , of 
F({n m , 0); y* ) = 

v v 

n ~k Jc m 1 — f n tn 

(yo YmO + n } + Xl 1 I Y ml * (A7) 

m=l m=l 

•k * 

where we have replaced P nr (y ) by F^n^.r^}; y ), y = Y m ^)> and shall replace 

n nr by n({n m , r m )). When this eigenfunction is substituted into the multimode 
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version of the integral equation (25), we obtain, much as before. 


(xj + yj” x 0 ) 


0 


P 3 * n m 
m - n m 0 ml 


n 


-A(xi y i + x 0 ) 


n 


(O'* * 

P ;(yo 3 m0 + yi } 


* n 
* m 


m" n m 


ml 


= n((n m , 0}) F((n m , 0); y ) , 


where 


n = V n 
/ ; m 


( 48 ) 


(49) 


m 


is the total number of photons in the states contributing to the given eigen- 
function. By equating coefficients of like terms on both sides of (48), we 
obtain a pair of iinear equations like those in (34), except for the replace- 
ments 


P n0 1> P<1)(( V ° }) 


= rL ( » .n 

V n m 0 V 


(i - v ln )(l - 


(0) p (0) ({n }) = PI P (0) = PI (1 - v )v_ n m , 

n v m m' 1 ' n A m 1 ' 0m 0m 

1 m 


v. = h.x N /(h.x N + 1) , j = 1,2, 

jm jms jms ’ J 

v_ = X N /(x N + 1) . 

0m m s m s 


(50) 


As the time-bandwidth product WT increases, the temporal eigenvalues X 


m 


become smaller and smaller, and we can approximate the denominators in the ex- 
pressions for v. , v, , and v_ by 1. Then for v >> WT >> 1, 

Um lm /m 


P (0) ((n }) - 
m 


P (1) ({n ro , 0}) 


I 1 exp(-x N) (xN) Ilm = exp(-N ) I I 
m ms ms s 

PI exp(-h lX N ) 


n_ 

exp(-N ) I I (x N ) , (51) 


m s 


exp(-h 2 X N ) (hjX N ) 

*- m c A m c 


■ h l" “P(- N s> U W”" ■ h l" p<0) «-.»• 


(52) 
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since hj + h 2 = 1 and the eigenvalues x m sum to 1. 

From (35) we see that the eigenvalue n({n , 0}) is now proportional to 

m 

P^^({n }) through a factor n n ' depending only on the sum n of the n m ' s > 

n({n 0}) = D ' p (0) ({n }) ; 

m n m 

that factor is a root of the equation 


n n 2 " (h l“ " A) V _ AhlI1(1 “ y l 2n) = ° * 


Furthermore, by the new version of (34) arising from (48), the coefficients 
xo and xi stand in a ratio depending only on the sum n and not on the individual 
n 's. If we go through the rest of the calculation of the error probabilities. 


m 


,(D 


we find that the sums in (39), (42), and (44) consist of terms with P v 1 ({n , 0)) 
or P^^({n }) multiplied by factors depending only on the sum n. We can there- 
fore combine terms having the same value of n = V n , replacing P^({n }) 

L m m m , 

when WT >> 1 

by the total probability of n photons under hypothesis H 0 , which/is the Poisson 

ability 
n D (0) 


probability P^^ in (45), and by virtue of (52) replacing the terms P^^({n , 0}) 


by hi n P 

We have used a scalar theory of the electromagnetic field, which seems to 
require that the light be linearly polarized. Unpolarized light can be divided 
into two statistically independent linearly polarized components, each of which 
is then broken up into temporal modes. The only change in our analysis is a 
doubling of the number of temporal modes, and as the number of such modes event- 
ually goes to infinity, our results for WT >> 1 must hold for unpolarized light 
as well. 
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The Error Probabilities 


We evaluate all the probabilities of error for observation over a rec- 
tangular aperture. In Fig. 1 is plotted, as a function of the parameter 

o = ca/AR and for various values of N , the average error probability for 

s 

choices between hypotheses Hq and , equal prior probabilities being assigned 

to each. As the two point sources separate, this error probability approaches 

\ exp(-N g ). When the sources are very far apart, the observer can count 

the numbers of photons in the planar mode no ' ( r ) and in all modes orthogonal 

to it. Finding any at all in the planar mode, he chooses H 0 ; finding any at 

all in the rest, he chooses Hi- The only possibility of an error arises when 

no photons at all are counted, whereupon he chooses Hq and Hj at random with 

equal probabilities 1/2. The probability of this event is expC-N^). The 

limiting value of % exp(-N ) is reached when a = 2, with only small deviations 

s 

above it for a > 2. If we say that when o = 2 the two point sources are resolved 
as well as they ever can be for a given average number N g of received photons, 
we require a separation of 2AR/a, which is twice that specified by the Rayleigh 
criterion. 

In Fig. 2 is plotted the detection probability versus the false-alarm 

probability Q 0 for N =2 and various values of a = ea/AR. These curves, along 

s 

each of which A is a parameter, are sometimes called the operating character- 
istics of the system. The portions between the dashed lines, where A = 1, are 
straight lines and represent the use of a randomized strategy each time no 
photons at all are observed, hypothesis Hj then being chosen with a certain 
probability f and Hq with probability 1 - f. As f varies from 0 to 1, the 
straight lines are traversed from left to right. To the right of the dashed 
lines is the region A < 1 (5 < 1/2); to the left is the region A > 1 (5 > 1/2). 
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Zero false-alarm probability can be attained with a finite detection 
probability 

Q d = 1 - exp[-N g (l - hxy! 2 )] = 

1 - exp{-N g (l - sine 2 ho)}; Qo = 0, A = ». (53) 

The strategy achieving this chooses hypothesis Hj if any photons are observed 
in any mode orthogonal to the planar mode rig'(r). Since this is impossible when 
a single source is present (Hq) , Qo =0. An error occurs under hypothesis Hi 
only when no photons happen to appear in those orthogonal modes. The average 
number of photons in the planar mode under Hj is 


E( n o|Hi) = 



cos y(x^ - x 2 ) d 2 rjd 2 r 2 = N h-^ 2 . (54) 


The average number in the rest of the modes is N g (l - hjyjZ)^ and the proba- 
bility that no photons are observed in them and H 0 is chosen is 
exp[-N g (l - t^yj 2 )], whence we obtain (53). The complementary probability 
1 - Q d has been plotted in Fig. 3 versus a for various mean numbers N g of 
received photons. A larger detection probability can be attained by accepting 
a positive false-alarm probability and adopting a much more complicated strategy. 

Infallible detection (Q d = 1) can be achieved by accepting a false-alarm 
probability 


Q o = exp(-N g y 0 2 ), Q d =1, A = 0, 

where yg 2 = 1 - yj? is given by (17). For a << 1, yg 2 = (Tio) t4 /720; thus the 
false-alarm probability is close to 1 unless N g is very large or the sources 
are well separated. To achieve this pair (Qq , Q d = 1) hypothesis Hi is chosen 
when there are no photons in the zero mode n 0 (r) given by (11); otherwise Hq 
is chosen. The mean number of photons in the zero mode under hypothesis H 0 is 
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given by (14) as N y 0 2 , and the probability that there is none is exp(-N y 0 2 ). 

s s 

Under Hj there will never be any photons in the zero mode, and it will always 
be correctly chosen. 


The Circular Aperture 

We list here for reference the formulas for the constants appearing in 
our solution when the aperture is circular, centered at the origin and having 
radius a. With m = pa = kea/2R = irea/AR, 

y 1 2 = 8m -1 [m + J x (2m) ] -1 [J^m)] 2 , 

h x = [1 + m" 1 J 1 (2m)], 

y 0 2 = 1 - y! 2 , h 2 = 1 - hj, 

h^! 2 = 4m _2 [J 1 (m)] 2 ; (55) 

the last describes the Airy pattern for a circular aperture. Here Ji(x) is 
the Bessel function of order 1. 
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3. One Point Source or the Other? 


An optical communication system could transmit binary information by 
turning on one or the other of two point sources for intervals of duration T. 

If the sources are very close, a distant observer will make errors in deciding 
which one is on during a given interval. Both sources can emit radiant power 
P s ,and they are separated by e. We can put one at the origin u = 0 without 
loss of generality when the sources are close to the optic axis normal to the 
aperture at its center (e << R) . Under hypothesis Hq the emitting source is 
at u = 0, and the spatial covariance of the aperture field is that given by 
(2); under Hj it is at u = e , and the spatial covariance function is 

< Ps 1) ^l> r 2 ) = 

(E S /AT) exp [i k c • (r! - r 2 ) /R] Frj Fr 2 * (56) 

Otherwise the sources are as before. The first task is to find suitable spatial 
modes for the aperture field. 

The planar mode n 0 '(r) = A 2 Fr is an eigenfunction of (pg^ (rj , r 2 ) ; an 
eigenfunction of <p g ^(ri, r 2 ) is similarly 
A exp(i k e* r/R) Fr . 


We therefore take as our basic modes n Q ' (r) and 
V(r) = [Ci' + Cj" exp(ip'x)] Fr, 
y' = ke/R, e = |e| , 


(57) 


choosing C \ ’ and C^" so that Hj ' (r) and n 0 ' (r) are orthogonal and properly 
normalized. We find 


Ci" = [(1 - q 0 2 )A)"" 2 , C] ' = -Ci" q 0 , 


qo = A -1 / exp(iy'x)d 2 r . 

>h 


(58) 
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The aperture A may have an arbitrary shape, and without loss of generality q Q 
can be taken as real. The remaining spatial modes are made orthogonal to 
^0 ' (r) and n i ' (r) by the Gram-Schmidt process. They are unexcited by the 

source under either hypothesis and can be disregarded. As before we start by 
assuming that WT << 1 and only a single temporal mode needs to be considered. 
Later we again extend our formulas to cover a great many statistically indepen 
dent temporal modes. 

The two spatio-temporal modes are represented quantum-mechanically by 
harmonic oscillators, and the density operators pq and under the two hypo- 
theses have Gaussian P-representations with mode correlation matrices given by 
(13) as 

( 0 ) „ . 

<P pr - N s 6 P 0 6 r0 » 

(!) 

<P pr = N s q p q r , 

(p, r) = (0, 1), qi = (1 - qo 2 )^ 2 ; (59) 

N g is the mean number of photons received from either source during (0, T) . 

The optimum strategy is again specified by the eigenvectors of pj - Apg 
which can be determined from an integral equation like (25) , except that only 
two modes instead of three are involved, and a single subscript suffices to 
label the eigenvalues and eigenvectors. Here 

R 0 (§ > y) = (1 - v 0 ) exp v 0 g 0 * Y o (60) 

“fc 

R i(£ » y) = (1 ~ v o) ex p t v o (90 + qi 3 i*)Cqovo + qiYi)]» 
v 0 = n s /(n s + 1) . (61) 

The eigenfunctions in the coherent-state representation of (24) are now 
* * n * * n 

F n Cy ) = x 0 Yo + x l Cqo YO + qi Yl ) , (62) 

which when substituted into (25) yields 
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( 63 ) 


p nt(xoqo n + x j) (qo^o + qi Bi ) 

-A(x 0 + XI q 0 n ) 3o n ] = n n p n (6 ) . 

P n = (1 - v 0 ) v 0 n , 

and equating coefficients of like terms gives the homogeneous equations 

p n qo 11 x 0 + ( p n " H n ) x x = 0, 

-(AP n + n n ) x 0 - AP n q 0 n x 2 = 0, (64) 

the determinant of which, when set equal to zero, provides the quadratic equa- 
tion for the eigenvalues, 

n n 2 - (1 - A) P n n n — A ( 1 — q 0 2n )P n 2 = 0. (65) 

When the prior probabilities of the two hypotheses are equal, A = 1, and 

n n = ± P n d " q 0 2n ) % - (66) 

The error probability, given by (21), is now 

oo 

p 0 + ^ V 1 - a - V n > !s ] 

n=l 



(67) 


An argument like that in Section 2 permits us to include multiple tem- 
poral modes in the limit WT >> 1 by simply replacing P^ by the Poisson probability 

P = N n exp (-N )/n! . (68) 

n s e s 

In Fig. 4 we have plotted the resulting error probability, postulating a rec- 
tangular aperture of width a in the direction parallel to the line between the 
sources. The separation e is embodied in the parameter a = ea/XR, and qg = 

sine a. The limiting value is again % exp(-N ), attained when a = 1, which 

s 

corresponds to the image separation prescribed by the Rayleigh criterion when 
a diffraction-limited optical system is used. Now, however, the image of only 
one point source will be present at a time. 
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Appendix 

Some Coherent-State Calculus 


Let a = (aQ } ai,...) be a column vector of coherent-state amplitudes, let 


ot + = (ag*,cti*, . . .) be its Hermitian conjugate row vector, and let K be an n x n 
Hermitian matrix. Then the multivariate Gaussian integral in its most con- 
venient form for calculations with the coherent-state representation is 


^ ... ^ exp(-a + Ka + 8 + a + a + y) j (d 2 ct-jL/ tt) = 


(det K) -1 exp(S + K -1 y) , 


(69) 


where d 2 a^ = da^ x da^, 8 and 1 are constant column vectors, and the integra- 
tion is carried out over the entire 2n-dimensional space of (a. , a. }. 

ix ly 

In order to derive the coherent-state representation of the density 
operator (18), as given in (26), one applies equation (9.11) of Glauber's 
paper to (18), with c = <p, 


R i(§> = l det ¥ 


' ^ . .^exp (~a +< £ *a + 8 + a + a + y - a + a) (d 2 ou/Tr) 


- |det (eg + I)j -1 exp[|3 + (cp 1 + I) 1 y] (70) 

by (69) with K = <p -1 + I, where I is the identity matrix. 

The integrals required when substituting the eigenfunctions of (32) 
into (25) are most easily derived from the generating function 

rn 


f i(8* , k) = 


f R i(§*> exp V 1 ^ K i Y i* “ i "Yi I | (d 2 y i /n) 

' U =0 -1 i =0 


= | det (cp (i) + I) 



i=0 

exp[0 + (cp^ + I) -1 


+ I + S " I + y] | | (d^i/r) 
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| de t (cp (i) + I) | -1 exp[6 + (cp -1 + I) 1 k], i = 0, 1. 

(71) 


We see that the effect of the integration is to replace y by k. Thus 
f l(§*> £) = (1 ~ vj ) (1 - v 2 ) expCvjBi^Ki + V2 62 * K 2 ) > 


(72) 


f 0 (§*. !5) = (1 - v 0 ) exp[v 0 (y 0 B 0 * + yi3i*)(yo<o + y 1 K 1 ) J - (73) 


Putting tco = 9 yo> K 1 = yi» k 2 = 0 in these, we get the generating function 

£t JJJ *“ ' *’ 

i = 0, 1, (74) 

f 1 f (§* * 9) = (1 - vj ) (1 - v 2 ) exp(cp viyi3i*) (75) 


Y .7 0 > ^1 — t 7 y 1 > k« 2 “ u J-ii Liic&e , we get uue geueidiiug luuuiuu 
-i* (g* » «P) = ^ ex P [cp(y 0^0* + y lYl*) - I Y | 2 ] | | (d 2 y i /iT) , 


with 


fo'(§*» V) = (1 “ v 0 ) exp[tp v 0 (y 0 Bo* + y 1 B 1 *) ] , 

since yg 2 + y^ 2 = 1. The required multivariate integrals are then (i 
• 2 
Y? n Ri(§*> y) ex P(~ 111 2 ) | | (d 2 Yi/ir) 


(76) 
0 , 1 ) 


i=0 


= 3 f.(§*, K)/3K! n | <=g , 


(77) 


(yoYO* + yiYl*) n Ri(§*, y) exp(-|y| 2 ) | | (d 2 y±/v 


) 


i=0 


- 3 n f i '(§*, c P)/ 9t P n | (p=0} (78) 


from which (33) follows immediately. Equation (63) is similarly derived. 
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Figure Captions 

Fig. 1. Average error probability P e in deciding whether two sources or one 

is present, versus the separation parameter o = ea/XR. The curves are indexed 
with the mean number N s of received photons. 

Fig. 2. Operating characteristics of the optimum strategy for deciding 

whether two sources or one is present. The curves are indexed with the separa- 
tion parameter o = ea/XR; N s = 2. 

Fig. 3. False-dismissal probability (1 - Q^) in decisions whether two 

sources are present, versus the separation parameter a = ea/XR. The false- 
alarm probability is zero. The curves are indexed with the mean number N g of 
received photons. 

Fig. A. Average error probability P e in deciding which of two sources 

separated by e is radiating, versus the separation parameter o = ea/XR. The 
curves are indexed with the mean number N s of received photons. 
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